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ABSTRACT 



Channel Estimation is an essential component in applications such as radar and 
data communication. In multi path time varying environments, it is necessary to 
estimate time-shifts, scale-shifts (the wideband equivalent of Doppler-shifts) , and the 
gains/phases of each of the multiple paths. With recent advances in sparse esti- 
mation (or "compressive sensing"), new estimation techniques have emerged which 
yield more accurate estimates of these channel parameters than traditional strategies. 
These estimation strategies, however, restrict potential estimates of time-shifts and 
scale-shifts to a finite set of values separated by a choice of grid spacing. A small 
grid spacing increases the number of potential estimates, thus lowering the quanti- 
zation error, but also increases complexity and estimation time. Conversely, a large 
grid spacing lowers the number of potential estimates, thus lowering the complexity 
and estimation time, but increases the quantization error. In this thesis, we derive 
an expression which relates the choice of grid spacing to the mean-squared quantiza- 
tion error. Furthermore, we consider the case when scale-shifts are approximated by 
Doppler-shifts, and derive a similar expression relating the choice of the grid spacing 
and the quantization error. Using insights gained from these expressions, we further 
explore the effects of the choice and grid spacing, and examine when a wideband 
model can be well approximated by a narrowband model. 
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CHAPTER 1 



Introduction 



Accurate channel estimation is an important component of digital communica- 
tion and radar. For our research, we concentrate on the estimation of multi path 
channels with sparse impulse responses. Such channels are encountered in many ap- 
plications such as high definition television (HDTV), communication near a hilly ter- 
rain, and underwater acoustic communication near the surf zone [1] [2]. Due to the 
sparse impulse responses of these channels, traditional estimation techniques such 
as least-squares result in over-parameterization and thus poor performance of the 
estimator [3]. Fortunately, the structure of these channels can be exploited using 
sparse reconstruction algorithms such as Matching Pursuit (MP) [4] [1], Orthogonal 
Matching Pursuit (OMP) [5] [6], Basis Pursuit (BP) [7], or Fast Bayesian Matching 
Pursuit (FBMP) [8] . Since these algorithms are better suited for the channels we con- 
sider, they give more accurate channel estimates than the traditional least squares 
method [3]. 

In this thesis, we narrow our focus to the case where the signals used for channel 
probing are wideband and the channel is rapidly time varying. In this case, the nar- 
rowband Doppler approximation used by other authors may not be valid [9] [2] . Our 
goal is to find the cost of using the narrowband Doppler approximation, as well gain 
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insight on the error resulting from this estimation strategy. In Chapter [2] we present 
a model for the channels we consider as well as a few channel examples. In Chapter [3] 
we present detailed explanation of the estimation strategy we use. In particular, we 
explain how in order to use a sparse reconstruction algorithm, a linearized model is 
constructed composed of either the effects of time-shifts and Doppler-shifts or the 
effects of time-shifts and scale-shifts (By scale-shifts, we mean the wideband equiva- 
lent of Doppler-shifts [9]). A design parameter of this linearized model is something 
we define as the grid spacing, which we describe in detail in Section 13.21 In general, 
the larger the grid spacing, the fewer the number of potential estimates the sparse 
algorithm has to choose from. In Chapter H] we derive an analytical expression which 
relates the estimation error to the choice of grid spacing when time-shifts and time- 
scales are used. Then in Chapter O we derive an analytical expression which relates 
the estimation error to the choice of grid spacing when time-shifts and Doppler-shifts 
are used. With these two expressions, we gain insight on the effects of the Doppler 
approximation as well as the effects of the choice of grid spacing. Finally, in Chap- 
ter [6] we present computer simulations of the estimation strategy. We focus on the 
application of digital communication, but we keep the problem general enough so 
that the results provide insight to radar systems as well. 
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CHAPTER 2 



Channel Model 

In general, the channels we consider are multi path channels with time varying 
delays and gains (see Figure [2]T] for a simple example). We define the channel impulse 
response as follows 

L 

Mt,r)4 J]/«(t)5(r-r«(t)), (2.1) 

where L is number of channel paths, T^'\t) is the delay in seconds of the l*^ path at 
time t, and f^'^t) is the complex gain of the l^^ path at time t. Over short periods of 
time, r'''(t) can be approximated by a first degree polynomial as follows 

T^'\t) = t''\0) + a^'H. (2.2) 

In (12. 2p . a*'' is the rate of change of the channel delay ^^-^f^, and r'''(0) is the delay 
when t = 0. By denoting c as the speed of propagation of the signal and 
length of the Z*^ path in meters, (12. 2p becomes 

c c c 

where C*'' is the rate at which the length of the l^^ path changes. 

One must be careful with joint estimation of the parameters r*''(0) and a'*'. As 
shown in [10], the parameters are coupled, and as a result, an estimation error in one 
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parameter may result in an estimation error in the other. This motivates us to use a 
different model for r^''(t) composed of uncoupled parameters. To do this, we will first 
note that at time t, the signal travels a distance of ct, while the path length x^'K^) 
is equal to + C^'^t- By setting these two distances equal and solving for t, we 

can denote 7*'^ as the time it takes the signal to travel from the transmitter to the 
receiver along the l^^ path 

ct = x"\0) + C'%=,ii) (2.4) 
^ (0=X^_ (2.5) 

Using the fact that C^'' = ca^'^ and x^^\^) = cr^'''(t), we can rewrite (12.51) as follows 

' c-ca« l-a«' ^ ■ ' 

Now by combining the results of (12.21) and (12. 6p . we can model the time varying delay 
in the following way 

r«(t) = (l-a«)7« + a'"t. (2.7) 

Since it has been shown that there is no coupling between the parameters 7''' and a'*' 
in [10], the goal of the estimator will therefore be to find a*'' (the rate of change of 
the delay) and 7*'' (the time it takes the signal to travel from the transmitter to the 
receiver). We can now rewrite the channel impulse response (12.11) as 

L 

h{t,T) = J] (7«(t)Vl - a«5(r - (1 - a^'^)Y'^ - a«t). (2.8) 

1=1 

By replacing with (7'''(t) vT^-a^ in (12. Sp . we in effect have normalized the model 
so that for any a^'\ the path is energy preserving when g^'\t) = 1 [9]. Now we will 
present two examples of time varying channels to further illustrate this model. 
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Figure 2.1: Multi path channel example 

2.1 Moving Receiver and Fixed Transmitter 

The first example is a single path channel with a fixed transmitter and a moving 
receiver (see Figure [2^2]) . In this example, the receiver is traveling with a vertical 
velocity of f{,'' and a horizontal velocity of v'f^\ The horizontal distance between the 
transmitter and the receiver is a^'^t), and the vertical distance is The path 

delay, T'''(t), can be expressed as follows 

,0... _ X'\t) _ Va(O^(t) + 0(O^(t) _ VKKO) + <t)^ + (0(0(0) +.;^'>t)^ 
^ c " c " c ^ ' 

We would like to approximate (12.91) to fit the model given in (12. 7p . To do this we 
will first take the derivative of (12. 9p with respect to t. 

ar«(t) (ff*"(0) + vllH)vll^ + (0("(o) + vi'Hyj^ 



dt 



;^(a«(0) +<t)2 + (0(0(0) +t;«t)2 



(2.10) 
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Now we will use the Taylor Series Approximation around t = 0. 



r'"(0) + — 



dt 

XW(0) a"HOK" + 0^"(O)^^" 
c cv/a(02(0) + 0(0^ 

c cx(')(0) 
X'" (0) v^l^ cos (0) + sin (0) 



+ 

c c 
Equation (12.141) follows from the facts that cos'(9^'^(t) 



X("(i) 



and sin'(9^'^(t) 



We can recognize that if we express a*'' and 7*'' as follows 

cos7?("(0) +^^*'' sin#"(0) 



,(0 
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(0 



x«(o) 



(2.11) 
(2.12) 

(2.13) 
(2.14) 

0(O(t) 

(2.15) 
(2.16) 



c(l - a«)' 

then (12.141) fits the model given in (12. 7p . with a*'' being approximately the rate the 
channel delay of the l^^ path changes and 7*'' being approximately how long it takes 
the message to travel from the transmitter to the receiver along the l^^ path. 



Receiver 



,(0 



,(0 



X'^KO),,^'- 



^(0(0) 



Source 



0("(O) 



Figure 2.2: Channel example (moving receiver and fixed transmitter) 
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2.2 Underwater Acoustic Channel 



Underwater acoustic channels near the surf zone provide a challenging application 
for sparse channel estimation techniques. Such channels are especially challenging 
due to rapid fluctuations in per-path gain and delay [11]. An example of a single 
path in such a channel is given in Figure 12.31 We define the tangent point (shown in 
Figure 12. 4p as the point on the water surface where the signal bounces on its path 
to the receiver. With the exception of the line of sight path, each of the other L — 1 
paths has its own unique tangent point. As the water surface moves, the tangent 
points move as well. We therefore expect each path delay to be time varying even if 
the transmitter and receiver are fixed. 

We now would like to find an expression for a single path delay r''^(t) that fits 
the model given in (12. 7p . To aid us in this task, we will use the model drawn in 
Figure 12. 3[ We assume that the tangent point moves with a horizontal velocity of 
vll^ and a vertical velocity of f^''. We can therefore express the horizontal distance 
between the transmitter and the tangent point as af{t) = a^g\0) + vl^H, and the 
vertical distance as (f)^'\t) = (p^'^O) + v!^H. Similarly, the horizontal distance between 
the tangent point and the receiver can be expressed as cr^'^it) = cr^'(O) — v^^H. We can 
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thus express the channel delay as follows 



rWU) = (2.17) 
c 



(2.18) 
(2.19) 



VK'(o) + <t)' + (</'^'Ko) + y^"t)' 



^/(<^<"(0)+^;i"t)2 + (a<J'(0)-<^)2 
+ -^^ . (2.20) 



In order to approximate fl2.20p to fit the model expressed in (12. 7p . we will use the 
Taylor Series Approximation as was done in Section [2TT1 

r'"(t) ^ r<"(0) + ^^'''^°^ t (2.21) 

ot 

x^'(o) , x'm 

c c 



a^'(O)< + 0'"(OK' . ^ 0'"(OK' + 4XOK . 

x!^^(o) + x^'(o) 



c 

^ cos^i;^(0)^;^'' + sin t^;^' (OX') + sin^^'(0)t;W + cos4'(0)t;fe ^ 23) 

c 

Equation fl^T^ follows from the facts that cos^9^'(t) = 4r^, sin^9^'(t) = 

sin'(9^'(t) = ^jTf^, and cos't?^'(t) = ^^y^. We now recognize that if a*'' and 7*'' are 



expressed as follows 



^(0 ^ cos<>(0)t;i'> + sin^g(0)^;<') + sin^^>(0)t;<') + cos<'(0)t;i'> ^^.24) 
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then (12.231) fits the model given in (12. 71) . From the figures of experimental data 
collected in [11], the linear approximation does appear to be valid over short periods 
of time. 




Figure 2.3: Acoustic channel model 
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Figure 2.4: Signal bouncing off the tangent point of the wave surface. 
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CHAPTER 3 



Estimation Strategy 

The objective of the channel estimator is to estimate the unknown channel pa- 
rameters {7^'', a'*', and accurately. We now present a general strategy to 
achieve this goal. The first step, which is outlined in Section [37T| is to send a training 
signal x{t) through the channel. In Section [3^ we describe in detail how the system 
model is linearized in the unknown parameters {7''' and a'''}^^. Then in Section [3l3] 
we demonstrate how a sparse reconstruction algorithm can use the received training 
signal and the linearized model to estimate the unknown channel parameters. For 
this chapter and Chapter IH we assume that the linearized model is composed of the 
effects of time-shifts and scale-shifts. Later in Chapter [5l we present modifications to 
this strategy to create a linearized model composed of the effects of time-shifts and 
Doppler-shifts instead. 

3.1 Training Sequence 

In order to estimate the channel, a sequence of pseudo-random QPSK pilot 
symbols is modulated, sent through the channel, and observed at the receiver. By 
pseudo-random, we mean that the values the pilot symbols are chosen to be random 
i.i.d., but they are known in advance at the receiver. The pilot symbols, which we 
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denote as {6„ G (±l±j)/-\/2}^^^^y25 transmitted using single carrier modulation. 
By using T to represent the symbol period, fc to represent the carrier frequency, and 
pit) to represent a square-root raised-cosine pulse, the passband transmitted signal 
becomes 



X 



N/2-1 

(t) = Re<je^W ^ bnp{t-nT)\. (3.1) 

n=-N/2 



The transmitted signal is then sent through the channel and observed at the receiver. 
The received passband signal r(t) can be expressed as follows 

r{t) = J x{T)h{t,t-T)dT + n{t) (3.2) 
= lx{r)J2fHmt-r-T^'\t))dr + n{t) (3.3) 



L 



= J2 ntMt - T^'\t)) + n{t) (3.4) 
1=1 

L ( N/2-1 1 

= P\t)Re \ eJ2-/c{t-(')W) ij^p(^t - r<')(t) -nT)\+ n{t), (3.5) 

1=1 [ n=-N/2 J 

where n{t) denotes white Gaussian noise. At the receiver, r{t) is low pass filtered and 
converted to baseband. We express the complex baseband received signal z{t) as 

z{t) = LPF {2r{t)e-^^^^^*} (3.6) 

L N/2-1 

= Kvit~nT -T^'\t))e'^^^f^^^'^^''^ +w{t), (3.7) 

1=1 n=-N/2 

where wii) denotes the noise after down conversion and filtering. Using equation 
(12. 7p . the complex baseband received signal z{t) becomes 

L N/2-1 

z{t) = J2 &„p((l-a''^)t-nT-(l-a''^)7«) 

1=1 n=-N/2 
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In practice, the complex baseband received signal is sampled at P times the symbol 
rate. Over the estimation period, G samples are yielded {za:}^=^~^ 

Zk = z{k^) (3.9) 

L N/2-1 

= Y^gl'Vl^^ Yl bnp{{l - a^'')k- - riT - {1 - a^'^W) 

1=1 n=-N/2 

where Wk = w{k^), g'j^'' = g'-'\k^), and H is the index of the first sample. 
We model g^^\t) using an A^^^ order basis expansion model so that 

9'\t) = Y,oTpm. (3.11) 

i=l 

or in the discrete case 

i=l 

Equation fl3.12p can also be written in vector form as follows 

5^'> = 5^0f'pi, (3.13) 

■t=i 

where gr^'^ and pi are G x 1 vectors. Using this basis expansion model, (13.101) becomes 

Nb L N/2-1 

Zk = Y.Y.^T Pi,k^^^' E W(l-«'')fcp-^T-(l-a"07'") 

i=l 1=1 n=-N/2 

^g-,2./.((l-a(0h(0W')fcf ) + (3,14) 

3.2 Linearization of the Model 

Even though they are uncoupled, estimation of the unknown channel parameters 
{7'''}^! c'-iid {a*''}^]^ is still challenging, because they are nonlinearly related to the 
complex baseband received signal samples {zk}'^^j^^^ ■ However, a linear approxima- 
tion can be made by modeling {zfc}^=^~^ as the summation of contributions from 
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many possible combinations of (7, a), but with only L of these contributions being 
"active" or non-zero. In particular, we consider values of 7 and a restricted to a 
uniform grid: 

7p = pA^ (3.15) 
aq = qAa (3.16) 

where p e {Nq, ...,No + - 1} and q G {-Na/2, Na/2 - 1}. The linear model is 
constructed of these N^Na possible combinations of 7^ and a^. Before describing the 
model further, it should be noted that A^^o, N-y, Na, A^, and A^ are chosen so that 

max 

2 ^ci — ^max 
(~^ l)Aa, ^ dmax 

where 7mm and 7max are the minimum and maximum values of 7 expected in the 
channel, respectively, and a^ax is the maximum value of a expected in the channel. 
Using the newly defined "grid" of restricted values ({7p}^^^^~"^ and {0'q}^2-Na/2) ' 
we can now write the linearized 03.14^ as follows 

No+Nr-l Na/2-1 Nt 

p=No q=-Na/2 i=l 

N/2-1 ^ ^ 

X XI bnP{{l - aq)k— - nT - (1 - ag)'jp) + ttgk—) + Wk - ek,{S.17) 

n=-N/2 

where is the error due to the linear approximation and Op g i is the channel gain 
corresponding to 7^, a^, and the basis vector pi. We now define as a N^NaNb x 1 
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vector composed of all the values of 9p^q^i concatenated together. Equation (13.171) can 
now be written as follows 

z = Ae + w-e. (3.18) 

In (13.181) . A is an appropriately defined G x N^NaNb matrix, where each column 
of A corresponds to a particular 7p, a^, and pi. The quantities z, w, e are G x 1 
vectors composed of the output samples, noise samples, and the linear approximation 
error samples, respectively. It should be noted that since the channel has L paths, 
we expect only LNh of the N^NaNb elements in 6 to be non-zero. Since 6 is mostly 
composed of zeros, we consider 6 "sparse" . In Section 13. 3^ we show how sparsity is 
used in the estimation of 0. 

3.3 Channel Parameter Estimation via Sparse Reconstruc- 
tion 

Our interest now turns to the sparse algorithm which is used to estimate the un- 
known channel parameters in (I3.18P using the complex baseband received signal 
z. In order to accomplish this estimation, the algorithm has complete knowledge, 
either implicitly or explicitly, of the G x N^NaNb matrix A, because A is composed 
of pilot symbols known at the receiver. A wide variety of sparse estimation algorithms 
can be used including Matching Pursuit [4], Orthogonal Matching Pursuit [5], Basis 
Pursuit [7], Greedy Basis Pursuit [12] and Fast Bayesian Matching Pursuit [8]. As 
explained in the introduction, these algorithms have an advantage of being more ac- 
curate than traditional estimation methods. This is because these algorithms exploit 
the fact that 6 is sparse (composed mostly of zeros). Additionally, these algorithms 
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have the advantage being able to work when the problem is underdetermined (i.e. 

G < N^NaNt). 

To allow further analysis of the estimation strategy, we now outline analytically 

how we expect the sparse estimator to behave. We assume that the channel gains are 

time invariant (i.e. A^^ = 1 yielding the simpler notation g^'^t) = ^*''). In this case, 

the estimator makes estimates of the channel parameters {7''\ a''\ ^^'^j^^. We assume 

that |6"''p = 1 so that the channel is energy preserving. We also assume that the 

number of significant paths the estimator finds L is equivalent to the true number of 

paths L. For each of the L estimated paths, parameter estimates are made which we 

denote as {7^'', a''\ and 0^''^)}f^i- Using these parameter estimates, a reconstructed 

complex baseband received training signal z{t) is made, which is expressed as 
L N/2-1 
z{t) = ^^^(''v^T^ J2 Kp{{i - a^'^)t - nT - {1 - a^'^)f'^) 

1=1 n=-N/2 

^Q-j2nM(l-amnW+aWt)_ ^3_^g^ 

We define the reconstruction error ez{t) as the difference between the actual com- 
plex baseband received training signal z{t) and the reconstructed complex baseband 
received training signal z{t): 

L 

e,{t) 4 z{t) - z{t) = J2 0^'^s^'\t) - ^«s«(i) - w{t), (3.20) 

1=1 

where s'''(t) and s^''(i) are defined as 

g(0(^-) A g-i27r/,((l-aW)-y(0+a(0t)y^-3^ 
N/2-1 

X Kp{{l - a^'^)t - nT - {1 - a^'^)Y'^) (3.21) 

n=-N/2 

A g-j27r/e((l-aW)7(')+aWt)^^ _ 
N/2-1 

X M(l-a''')*-^7^-(l-«''')f'')- (3-22) 

n=-N/2 
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The quantity ||e^(t)||^ can be upper bounded as follows 



|e.(t)|r 



1=1 



< ^||^<"s<"(t) -e'''s«(t) +\\w{t)\\ 



(3.23) 
(3.24) 



where the operator 1 1 1 1 is defined as follows 



\a{t)\f ^ I \a{t)f dt. 



(3.25) 



We assume that the sparse estimation algorithm selects channel estimates 
{7''', a*'', ^*''}^]^ that approximately minimize ||e2(t)||^. In particular, we assume that 
for each channel path, 7*'' and a*'^ are selected as the values closest to 7''' and a*'', 
respectively, on the uniform grid fl3.15p - fl3.16p . Therefore, we can say that a''^ and 7''' 
are at most half a grid space away from a''^ and 7''^ respectively, or in other words. 



I 1-2 



(3.26) 
(3.27) 



For each channel path, the component s'''(t) is yielded from 7''' and a*''. Given this 
s^'^(t), we assume that the sparse estimation algorithm minimizes the upper bound 
of (13.241) by choosing each value of as follows 



e 



(0 



(g'"(t),sW(t))gw 
(sw(t),sw(t)) ' 



(3.28) 



where the operator ( , ) denotes the complex inner product 



(a(t),6(t)) 4 / a{t)h*{t)dt. 



(3.29) 
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Equation fl3.28p can be simplified further by noticing tliat (s'''(t), s'''(t)) reduces as 
follows 

{s''\t),s''\t)) = / v^r^a^ Yl M((l-a*'')^-rT- (l-a''))f' 

r=-N/2 

N/2-1 

xVl J2 blp{{l - d'^'^)t - qT - {1 - d^'^)f^) dt 

q=~N/2 

= 5Z E Kb; piv)piv + rT-qT)dv (3.30) 

r=-N/2q=-N/2 
N/2-1 N/2-1 

= E E M;i?((g-r)r) (3.31) 

r=-N/2 q=-N/2 
N/2-1 N/2-1 

= E E ^'■^^^-'^ (3.32) 

r=-N/2 q=-N/2 

= N, (3.33) 

where -R(t) is the raised-cosine impulse response and 5„ is the Kronecker delta func- 
tion. Equation 03.311) follows from the fact that two square-root raised-cosine pulses 
integrate to raised-cosine pulse as follows / p{t)p{t — t) dt = R{t). Equation 03.321) 
follows from the fact that R{t) is equal to when r is a nonzero multiple of T, and 
equal to 1 when r = 0. Finally, 03.33P follows from the fact brb^ = 1 if r = g. Plugging 
in the results of 03.331) into 03.281) . we find that 

^,0 ^ {s^^\t),s^^{tW^ _ 
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CHAPTER 4 



Estimation Error Using a Time-Shift/Scale-Shift Model 



4.1 Error Metric 

To analyze the performance of the estimator, we need to measure the estimation 
error. There are many ways to do this, one of which is to directly compare the true 
channel parameters {T^'\a^''\0^'''}^^-^ with the parameter estimates {f*'', a*'', 
However, since it is unclear how to assign the relative importance of these parameter 
estimation errors, this may not be the best way to proceed. Instead, we measure the 
estimation error by evaluating how well on average the channel parameter estimates 
allow us to reconstruct a random data signal. To do this, we consider a series of M 
QPSK random i.i.d. data symbols {cn}^£^j^/2 modulated the the same way as the 
series of N training symbols {bn}n=-N/2- ^ passband signal modulated from these 
data symbols travels through the same channel as the training symbols so that the 
actual baseband received data signal y{t) can be written as 

L M/2-1 

y(t) = ^^('^v^r^^ ^ c^p((l - a('))i - riT - (1 - a^'OT^'O 

/=1 m=-M/2 
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Using the parameter estimates, the reconstructed baseband received signal y{t) can 
be expressed as 

L M/2-1 

y(t) = ^^^v^T^ CmPiil - a^'^)t - nT - {1 - a^'^m 

1=1 m=-M/2 

X g-i27r/,{(l-aW)7(0+a(0t)^ ^^^2) 

We define ey{t) as the difference between the baseband received data signal and the 
reconstructed baseband received data signal 

L 

ey{t) 4 y{t) - y{t) = - (4.3) 

1=1 

where 

A e-.2^/c((i-a(0h(0W')t)^/i3^ 

M/2-1 

X CmP{{l - a^'')t - riT - {1 - a^'^)Y'') (4.4) 

m=-AI/2 
M/2-1 

X 5^ c^p((l-a«)t-nT-(l-a«)7^'')- (4-5) 

m=-M/2 

Since we are interested in how well the parameter estimates allow us to reconstruct 
an unknown data signal, we use -^E {||ej/(t)||^} as the error metric, where E{ } 
denotes the expected value. The reason we divide by M is because we expect | |ej^(t) 1 1^ 
to increase as M increases. Thus, -^E {\\ey{t)\f^ can be thought of as a metric of the 
estimation error per data symbol. Generally the larger the value of -gE |||ey(t)||^}, 
the worse the estimator performed. In the following sections, we find a relationship 
between the error metric -^E 1 1 |ej,(t) 1 1^} and the grid spacing (A^ and A^) defined 
in Section 13.21 We first examine the case when the channel is time invariant in 
Section 14.21 while in Section 14.31 we consider the more complicated time varying 
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case. Then, in Section 15. 2[ we consider the resulting error when the narrowband 
approximation is used. 

4.2 Average Estimation Error of LTI Channels 

The objective of this sections is to find a relationship between the grid spacing 
(Aq and A^) described in Section [3l2] and the estimation error metric -^E {||ej^(^)||^} 
described in Section 14. 1[ Considering a LTI channel (as opposed to a time varying 
channel) greatly simplifies this problem. In particular, the channel impulse response 
becomes 

L 

/i(r) =^^(')5(r-7'')), (4.6) 
1=1 

and the complex baseband received training signal becomes 



N/2-1 



:{t) = J2 e''^e~^^^f^'<^'^ h^p{t -nT~ ^''^) + w{t). 

n=~N/2 



(4.7) 



1=1 



Additionally, in the time invariant case, (14.41) and (14. 5 p can be simplified as follows 

M/2-1 



'j2nfc^,W 



J2 c^pit - riT - Y'^) 



-j27rf^jW 



m=-M/2 
M/2-1 



J2 CmP{t-nT -Y'^). 



(4.8) 
(4.9) 



m=-M/2 

From (14. 3p . if follows that the error metric can be upper bounded as follows 

2 ■ 



1=1 



< 



(4.10) 
(4.11) 



1=1 

\2 



Instead of solving for -^E {\\ey{t)\\ } in terms of A^ directly, we instead solve for 
the upper bound in (14. lip , thus finding an upper bound on the error metric. Focusing 
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on a particular value of / G {1, 2, L}, we find (as shown in Appendix 

2 



M 



-E 



2|^(" 

w 

M 



-E 



^(0 



dt 



Rei—S'\t)d^'>{t)} dt 



\d''\t)\^ dt 



(4.12) 



We now evaluate each of the three terms in (14.121) separately, and then combine the 
results. 



The first term of (14. 12p becomes 

2 



-E 



5(0 



-E 



^(0 



0(1) 



E 



dt 



S'\t) 



dt 



(4.13) 



|^<"lMl(.«(t),S<"(t))|' 



M 



-E 



M/2-1 



M/2-1 



xeJ /" J2 c,p(t - rT - 7<") J2 c*gp{t- qT -j'^'^)dt\ . (4.14) 

[^"' r=-M/2 g=-M/2 J 

In (14.131] , we were able to split the expectation into the product of two expectations, 
because the values of the pilot symbols {&n}^=l~/v/2 independent of the values of 
the data symbols {cm}^=lM/2- simplify the second expectation in (14.141) . we note 
that, as was explained in Section 13. 3^ the integral of two square- root raised-cosine 
pulse can be reduced as follows 



j p{t -rT- f'^)p{t -qT - 7'") dt = R{{q - r)T) = 5r-q. 



(4.15) 
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Using (14.151) . we can simplify the following expression 

A//2-1 M/2-1 



I Crp{t-rT-r') clp{t - qT - Y') dt 

r=-M/2 r=-M/2 
M/2-1 A//2-1 

= ^ ^ Crcl p{t - rT - Y''>)p{t - qT - f'^) dt (4.16) 

r=-M/2 r=-M/2 
M/2-1 A//2-1 

= E E ''rCldr-, (4.17) 
r=-M/2 r=-M/2 

= M, (4.18) 

where (I4.18P follows from the fact that CrC* = 1 if r = g. From (I4.18p . it follows that 
the second expectation of (14.141) is equal to M. Before reducing the first expectation 
of firaj) . we note that s<'^(t) and s^'\t) (defined in flX^ and ([3221), respectively) 
can be reduced as follows for the case of a LTI channel 

N/2-1 

s«(t) = e-^'2'^^^^'" Y bnP{t - riT - Y'^) (4.19) 

n=-N/2 
N/2-1 

s«(t) = e-^''"^^^''' Y bnP{t-nT-f'^). (4.20) 

n=-N/2 

Keeping this in mind, we express the first expectation of (14.141) as 

|2 



E 



— ^ E )^Khlhlh^ j pit -tT- Y')p{t -qT- f 'y) dt 



p{t-nT - -f^'^)p{t -niT - 7«) dtj | (4.21) 



q,r,m,n=—N /2 

X 

^ N/2-1 
q,r,m,n=—N /2 

E {hrhXh^R{{q - r)T + 7''^ - 7'")^((m - n)T + 7''' - 7''')} • (4-22) 

To simplify (14.221) further, the expected value of the product of four random variables 
brb*b^bm has to be evaluated. To do this, we note that because the pilot symbols 
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{bn}n=^N/2 pseudo-random QPSK i.i.d. symbols, E{6„} = 0, E{6„6m} = 0, and 



r = m 



E{^n^m} = ^n-m- It therefore foUows that E [brb*b'^bm} is equal to 1 if g 
n, if q = r ^ m = n, OT if q = m ^ r = n. In all other cases E {^g^r^m^n} equal 
to zero. Also notice that in f l4.22p . the case q = r = m = n occurs times, the case 
q = r ^ m = n occurs A^(A^ — 1) times, and the case q = m ^ r = n occurs A^(A^ — 1) 
times. Therefore, f l4.22p can be reduced as follows 

N/2-1 

J2 ^ {hrb%bmR{{q - r)T + 7<" - 7'")i?((m - n)T + 7^'' - 7* 



q,r,m,n=—N j2 

- ^{A^i?'(7^" - 7'") + N{N - 1)R\Y'' - 7'")} + 



(0 



(4.23) 
(4.24) 
(4.25) 



where 



N/2-1 N/2-1 



a 



(0 A 



(4.26) 

q=-N/2 q=ir=-N/2 

By combining the results of fl4.18p and fl4.25p , we find that first term of (14.120 can be 
expressed as 



1^ 



(01 



M 



-E 



0W 
0W 



dt 



(4.27) 
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We now simplify the second term of fl4.12p using similar arguments 

2 1^^'' 



M 



eJ /Re 
V 


(eh 


M 1 




MN 


E{( 


MN 1 


H 


C A//2- 


-1 


xE E 


q,r=-M/2 



(4.28) 

E<| / S'\t)d^^>{t)dt\ \ (4.29) 

(4.30) 



N/2-1 

J2 b^b:R{{m-n)T + Y''-¥"] 

m,n=-N/2 



_ (4.31) 

,r=-M/2 J J 

= ' i?(7"' - f '')^(7''' - f (4.32) 

= 2|e«|'i?2(yo _^(0) (4_33) 

Equation fl4.32p follows from the facts that EjcgC*} = 6q-r, E{bmbn} = Sm-n, and 
R{t) is real for all r. 

Using the arguments we used to reduce the first two terms fl4.12p . the third term 
of (14.120 reduces as follows 

(4.34) 

Finally, combining the results of K27} . K33tl . and (14:31) into fl412|) . we find that 
1 



-E 

M 



In Appendix [B| it is shown that by assuming — 7*''| < ^, a''^ can be upper 
bounded as follows, a''' < Using this fact and keeping in mind the channel is 



25 



energy preserving Yll'=i = 1, we can combine fl4.35p with (14. lip to find 

^^{\\ey{t)\n < t.\^'f(l~R^l'^-r) + ^-^) (4.36) 

1=1 ^ ^ 

n 6Q4 

= l-i?2(yo_^(0) + l!:^. (4.37) 

Recalling our earlier assumptions that — 7*'^| < ^ (assumed in Section [3l3|) and 
1^(0 — ^c)] < ^ (assumed in Appendix [B|) . it follows that must be selected so that 
< -J. Keeping in mind that R{x) = -R(|x|) and -R(|a;|) decreases as \x\ increases 
for |x| < f , fOTI) can be upper bounded as follows 

^E{||e,(t)||^} < l-i?%«-f'>) + ^ (4.38) 

5,A^, 0.694 
< 1 - i?2(_I) + __. (4.39) 

This is the expression we were looking for to relate the choice of grid spacing A^ with 
the estimation error metric {||ey(t)||^}. Fig. 14. II shows a plot of 1 — i?^(^) + ^^^ 
vs. A^ on a T normalized scale, when N = 100. As can be seen, the finer the grid 
spacing, the lower the average estimation error. 

4.2.1 Simulation of LTI Channel Estimation Error 

To verify the results in (I4.39p . we ran a simulation in Matlab of a one path LTI 
channel using a carrier frequency of 18 kHz, a symbol rate of 24 kHz, and a sampling 
rate of 72 kHz (3 times the symbol rate). The channel had unity gain and a time 
invariant delay of 7 ms. The simulation followed what was outlined in Chapter [31 
Specifically, a signal composed of a series of 100 training symbols (A^ = 100) was sent 
through the channel and an estimate of 6"^-' was made based on (I3.34p . After that 
a signal composed of 100 data symbols (M = 100) was sent through the channel. 
The error per data symbol ||ey(t)||^ was calculated for 4 different estimates. The 
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Figure 4.1: Plot of 1 - R^{^) + ^ vs ^ 

simulation was repeated 30 times and the average results are shown in Fig. 14. 2[ As 
expected, the estimation error increases when the the grid spacing increases. It 
can also be seen that the upper bound of f l4.39p was satisfied for all values of 

4.3 Average Estimation Error of Time Varying Channels 

In Section 14.21 we found a relationship between the estimation error metric 
-^E {I |ey(t)| 1^} and the grid spacing A^ for the case of linear time invariant channels. 
We now consider channels with path delays that vary with time (i.e. a^'' may not be 
equal to zero for / G {1,2, ...,L} ). We still assume that the channel path gains are 
time invariant (i.e. g^'^t) = 6"'' for / = 1,2, ...L), the channel is energy preserving 
J2f=i 1^'*' P = 1) and that the number of estimated paths L is equal to the true number 
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Figure 4.2: Derived upper bound and average simulated error for the LTI case 



of paths L. Additionally, we now assume that the number pilot symbols is equal 
to the number of data symbols M. 

First we recall, as noted in Section that the error metric can be upper bounded 
as follows 

>{im*)ip}<^i:e{ 

As we did in Section 14.21 we evaluate this upper bound of instead of evaluating 
^E{||ey(t)||^} directly. We again focus on a particular value of / G {1, 2, L} and 



(4.40) 
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use the fact that (as shown in Appendix \M 



M 



(01 



M 



-E 



0(1) 



eh 



dt 



(01 



M 



-Ejy" \S'\t)\^ dt 



(4.41) 



Now we evaluate the 3 terms in (14.411) separately and later combine the results. 

Using (13.341) and the assumption = M, the first term of (14.411) can be expressed 
as follows 

2 



1^ 



(01 



M 



-E 



0ii) 



dt 



\0 



(01 



M 



-E 



^(0 



0(1) 



e\ / d''\t)d^^>{t)dt 



(4.42) 



1^ 



(01 



M 



-E 



(^'"(t),g">(t))) i [ (s'"(t),s''>(^)) r 



M 



I 



M 



M/2-1 

/ v^r^Iw Yl ^-^^^^ ~ ^''')^ - rT - (1 - a«)7'") 

r=-M/2 



M/2~l 



X 



v^r^a^ 5Z c;p((l-a«)t-gT-(l-a«)7''')^^i 

g=-M/2 



(4.43) 
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To simplify the second expectation of (14.431) . we note that by making the substitution 

V = {1 — d*'')t — rT — (1 — a^'^)7*'\ we can reduce the following expression 

M/2-1 

VT^^ J2 c,p((l-a''')t-rT-(l-a«)7'") 

r=-M/2 
M/2-1 

xVl -a(') c^piil - a''^)t - qT - {1 - d^'^)f'^) dt 

q=-M/2 

= r^l^ 5Z 5Z Crc; p{v)p{v + rT-qT)dv (4.44) 

r=-M/2q=-M/2 
M/2-1 M/2-1 

= E E ^^'^^-^ (4.45) 

r=-M/2 q=-M/2 

= M. (4.46) 

It therefore follows that the second expectation of (14.431) is equal to M. Equation 
(I4.45P follows from (14.151) . and (14.461) follows from the fact that c^c* = 1 if g = r. 
Before simplifying the first expectation of (14.431) . we use the fact that (s^''(t), s'''(t)), 
as shown in Appendix [Cl can be approximated as follows 

M/2-1 M/2-1 

(s«(t),s"xt)) = E E 6,6;F^';^e^-2-^^«i-'^^'^)^^'^-(i-'^^'^)-'^'^) 

r=-M/2q=-M/2 

xW{qT -^ + {1- a'")(7<'' - 7'"), Z^'"), (4-47) 

where iy( ) is the wideband ambiguity function defined in (1C.7I) . F^^'^l is defined in 
(1C.4I) . and = Iz^w ■ Using (14.471) . we can express the first expectation of (14.430 
as follows 



E 



(sW(t),g'')(t))) ) / (s")(t),gW(t)) 
M J \ M 



M/2-1 



M/2-1 

E 

r,q,m,n=—M/2 
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where 



A« 4gT-^ + (l-a<'))(7")-7<")- 



(4.49) 



As we did in Section I4.2[ we use the fact that E \^brb*b^bm} is equal to 1 if g = r = 
n = m, ii q = r ^ n = m, or ii q = m ^ r = n. In all other cases E {^brb^b'^bm} is 
equal to zero. We also use the fact that F^llF!n]n = 1 if q = m and r = n to simplify 
(OH]) as follows 

( M/2-1 

\ brbib:b^F!^;iFji:>;w{xz, nwCK]m, n 



where 



r,q,m,n=—M/2 



(0 



( M/2-1 

\^r=-M/2 

M/2-1 M/2-1 

m=-M/2 nj^m=-M/2 

{M/2-1 M/2-1 

r=-M/2 q=-M/2 
M/2-1 

m,n=-M/2 
( M/2-1 

I / J m,m n,n " \^m,mi 1^ J" \^n,ni 

m,n=-M/2 

2 

^ Ml2-^ 



{^)^ 



UJ 



(I) 



M2 



M/2-1 
n=-M/2 



(0 



Af/2-1 



(0 A 



e^r=-M/2 



(4.50) 
(4.51) 



(4.52) 



(4.53) 



(4.54) 
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Now by combining the results of fl4.46p and (14.531) . we can express the first term of 
fOTjl as follows 



1^ 



(01 



M 



-E 



em 



0W 



dt 



N/2-1 
n=-N/2 



+ 



M 



1^ 



(01 



M2 



M/2-1 
n=-M/2 



(4.55) 



(4.56) 



Using some of the arguments we used to reduce the first term of (I4.4ip . the second 
term of (I4.4ip can be simplified as follows 



2\e 



(01 



M 



E|y"Re|^rf'''(t)rf«*(t)| dt 
-Re|E||^rf"'')(t)d«*(t)rft 



2\e 



2|gW|' 
M 



2\9 



(01 



M 



-Re<E 



(g'"(t),sW(t)) 
M 



E<^ / d^'\t)d^^>(t)dt 



2\9 



(01 



M2 



2\9 



(01 



M2 



2\e 



(01 



M2 



{M/2-1 M/2-1 
n=-M/2 m=-M/2 

{2 
M/2-1 
n=-M/2 

2 

M/2-1 
n=-Af/2 



(0^ 



(0^ 



?(0^ 



(4.57) 
(4.58) 
(4.59) 
(4.60) 

(4.61) 

(4.62) 



Equation (Ml follows from fOTll . (109ll . and the facts that E {bnb*^} = 5n- 

and E{CnC*^} = 6n-m- 
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Finally the third expectation of (14.411) can be reduced as follows 



E<^ / \S'\t)\^ dt 



M 



(4.63) 



Thus, combining the results of (I4.56p . (I4.62p . and (I4.63p . we can express (14.4 ip as 
follows 



M 



- 



3(01 



M2 



M/2-1 
n=-A//2 



(4.64) 



Furthermore, if we merge the results of (lOm and flTOIl . j^E {\\ey{t)\f} can be 
upper bounded as 



M 



E{||e,(t)ir} 



< 



El" 

1=1 



"M2" 



M/2-1 
n=-M/2 



(0^ 



+ |0"fc2;«. (4.65) 



Finally, by applying the results of Appendix [D] and Appendix [El (14.650 becomes 



(4.66) 



«=1 

X I sine 



+ 



2M 



1 - a(') 



47r 



2M 



1 - 



An 



M 



(4.67) 



where i? < vr is chosen to make sinc(x) > cos(i?x) for |a;| < |. Although this is a 
rather complicated expression, we can use it to gain insight on the performance of 
the estimator. For instance, if we assume that M > 100 and that a''' = .005 for all 
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I e {1, L}, and plug in a value of B = 1.8138, we find 

^E{||e,(t)||^} (4.68) 
< 1 - 0.2450 cos2(M160A,) 

X sine (.2488MAa (271 f^T + .5774)) 

\ 2 1 2011 

+ sine (.2488MA„ {^TrfcT - .5774)) j + -j^, (4.69) 

which follows from the assumption = 1. It can now be seen that 

|||ey(t)||^} decreases for decreasing and A^. Additionally, as the fractional 
bandwidth decreases, ^^gE {||ey(t)||^} increases. Next we examine how tight the 
upper bound of (14.671) is through simulation. 

4.3.1 Simulation of Time Varying Channel Estimation Error 

To examine the tightness of the bound in (14.691) . we ran a simulation of a one 
path channel with parameters a'^^ = —0.001, and 7*^' = 0. We sent a training signal 
modulated of 300 random symbols through the channel. This training signal had a 
fractional bandwidth = 1. Parameter estimates were given by 7*^' = 7'^' + ^ 
and a*^^ = a*^^ + Thus, we simulated the worst case scenario where the param- 
eter estimates were half a grid spacing away from the true parameter values. The 
estimate of was made based on (I3.34p . Next a data signal modulated from a 
series of 300 random data symbols was sent through the channel. The error metric 
-pE |||ey(t)||^} was calculated for different many choices of grid spacing. This simu- 
lation was repeated 10 times and the results were averaged. We set B = 1.8138 and 
plugged in the simulations parameters into (I4.69P to calculate the error upper bound 
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as 

^E{|Mt)ir} 

< 1.004- 0.2483 cos^( °-^°g'^^ ) (sine (193.32AJ + sine (106■8AJ)^ (4.70) 

Figure 14.41 plots this upper bound while Figure 14.31 plots the average simulated 
error. Figure 14.51 shows the two together, and it can be seen that the simulated error 
is less than the upper bound we derived. These results are also shown in Figure 14.61 
and Figure 14. 7[ Figure 14.61 plots the simulated error and the derived bound for fixed 
values of Aq, while Figure 14.71 plots the simulated error and the derived bound for 
fixed values of A^. It all cases it can be seen that the bound derived in fl4.67p does 
indeed hold and appears to be relatively tight. 




Figure 4.3: Average simulated error using a time-shift/scale-shift model 
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Figure 4.4: Derived upper bound for time-shift/scale-shift model 
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Figure 4.5: Derived upper bound vs. average simulated error for time-shift/scale-shift 
model 




Figure 4.6: Simulated error and derived upper bound vs. ^ for fixed values of 
using time-shift/scale-shift model. Simulated error is plotted with solid lines while 
the derived bound is plotted with dotted lines. 
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Figure 4.7: Simulated error and derived upper bound vs. for fixed values of ^ 
using time-shift/scale-shift model. Simulated error is plotted with solid lines while 
the derived bound is plotted with dotted lines. 
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CHAPTER 5 



Estimation Error Using a Time-Shift/Doppler-Shift Model 

5.1 B ackgr ound 

We now consider the case when a simphfying narrowband Doppler approximation 
is used to construct the hnearized model described in Section 13.21 Using this approxi- 
mation, the effects of time scahng by the quantity 1— a are approximated by a Doppler- 
shift by the quantity — 27r/ca as described in [9]. For example, consider a transmitted 
passband signal x{t) with a carrier frequency fc- If x{t) travels along an energy pre- 
serving single path channel with a time varying delay r(t) = {l—a)'~f+at, the resulting 
received signal r(t) = y/1 — ax({l — a)(t — 7)) is approximated as e~^'^-^''"''*~'^-'a;(t — 7). 

As we did in Chapter HI we assume that the path gains are time invariant 

= VT^-aJ^O'^'^''}^^^), the number of estimated paths is equal to the true 
number of paths {L = L), and the number of pilot symbols is equal to the number of 
data symbols {N = M). In the previous sections, we said that the sparse estimator 
makes estimates {a'''}^]^ of the true channel values where a'*' is the rate of 

change of the channel path. In this chapter, we say the sparse estimator instead 
makes Doppler-shift estimates denoted as {/'''}^^. We continue to use the tilde 
symbol to denote the Doppler approximated estimates. Thus, the sparse estimator 
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makes parameter estimates {7*'', 9^'^^}^^^, where 7'^'' and are the estimates of 7^'^ 
and 6^'\ respectively. Equivalently, we can say the sparse estimator makes parameter 
estimates {7''', where a*'' = y^. As outhned in Section [3721 a hnearized 

model is constructed and parameter estimates {7*'\S'''}^^ are grid restricted with 
grid spacing and A^. Using these parameter estimates, the pass band signal is 
approximated as follows. 



(t) = J2 ^*Vl - a('H{{l - a'")(t - 7'")) + nit) (5.1) 
1=1 

L 



1=1 

L 

1=1 

L [ N/2-1 



1=1 n=-N/2 

+ n{t) (5.3) 

We express the Doppler approximated estimate of the complex baseband received 
signal z{t) as 

z{t) = LPF{2f(t)e"^'2"^=*} (5.4) 

L N/2-1 

= ^^~(0e-^-2-/c({i-«)7(')+a(0i) ^ i,^p^t_^T-Y'')+wit) (5.5) 

1=1 n=-N/2 
L 

= ^^'')5«(t) +w;(t), (5.6) 

where 

N/2-1 

5(0(^) Ag-.2./.((i-a(0)^(0+a(0t) j2 b,,p{t-nT-f'^). (5.7) 

n=-N/2 
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For each estimate pair 7*'' and a*'', s*''(t) is yielded. Using similar arguments which 
lead to (I3.34p of Section I3.3|, we assume that the estimate of 6'*'^ can be expressed as 



We now need a way to measure the estimation error resulting from the use of the 
Doppler approximation with channel parameter estimates {7^'', a''', In Sec- 

tion |1]T], we considered a complex baseband received signal y{t) and another complex 
baseband received signal y{t) reconstructed from the channel parameter estimates. 
Both y{t) and y{t) were modulated from the same series of M PSK symbols. We 
defined ey{t) = y(t) — y{t) and used -gE |||ej,(t)||^} to measure the estimation er- 
ror. For this chapter, we define y{t) as the Doppler approximation estimate of y{t) 
constructed from the channel parameter estimates {7''\ o-'*', 6''*'};^;^ , so that 

L 

y{t) = Y.~^''d^'\t), (5.9) 
1=1 

where 

N/2-1 

Ag-J-2./c({l~a(0)^(0+a(0t) ^ 5^p(^_^2^_^(0). (5.10) 

n=-N/2 

We define ey{t) = y(t) — y(t) and use -pE |||ey(t)||^} as the metric of estimation 



error. In Section 15.21 we derive an expression which relates this metric to the grid 
spacing A-y and A^. Then in Section [573] we evaluate how this expression compares 
to the expression derived in Section 14.31 where the Doppler approximation was not 
used. This analysis helps us gain insight on the effects of the grid spacing choice and 
the effects of the Doppler approximation. 
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5.2 Average Estimation Error 



As we did in Section Upland Section IT3l we evaluate the upper bound of the error 
metric instead of the metric directly. 

L 



1=1 ^ 



(5.11) 



We again focus on a particular value oil G {1, 2, L} and find, by a similar derivation 
to the one presented in Appendix [Aj that 

2 



M 



-E 



-E 



0(1) 



0(1) 



dt 



\0 



(01 



M 



fE|y'Re|^J<''(t)rf«*(t)| dt 
-Ejy \S'\t)\^ dt 



(5.12) 



The task now is to evaluate the 3 terms in 05.121) and then combine the results. 
The first term of fl5.12p can be reduced as follows 

2 



-E 



0(1) 



0(1) 



S'\t) 



dt 



-E 



0(1) 



E^ / S'\t)d^'>{t)dt 



(5.13) 



f {s''Kt),PKm \ f (.(')(t),5'')(t)) 



M 
xE 



M 



M/2-1 



M 

M/2-1 



f J2 c,j9(t - rT - f ") J2 c^pit- qT -Y'^)dt\ . (5.14) 

r=-M/2 q=-M/2 J 

Using fl4.15l) . it can be shown that the second expectation in 05.141) reduces to M. 
To simplify the first expectation in 05.141) . we use the fact that (s*''(t), s*''(t)) can be 
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approximated as (shown in Appendix [F] ) 

A//2-1 M/2-1 
r=-M/2 q=-M/2 

xl^(gT-^j + 7«-7«,/3'''), (5.15) 

where Fgr is defined in flF.2p . and = 1 — a*'^ Using (15.151) . we can express the 
first expectation of fl5.14p as follows 



E 



M J \ M 



M2 

where 



1 I I 

-E ^ 5^ &.&;&;&..4.^;,„W^(A« , /3«)W^(A^:|„, ^ , (5.16) 

r,q,m,n=—M/2 



M = gT-;^ + f''-7'". (5.17) 



As we did in Section [4.21 and Section [4.31 we use the fact that E|6r6*6*6m} is 
equal to 1 if g = r = n = m, if g = r 7^ n = m, or if g = m 7^ r = n. In all other 
cases E {6^6*6* 6m} is equal to zero. Additionally, we use the fact that Fq^rF^^n = 1 
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a q = m and r = ?7, to simplify fl5.16p as follows 

M/2-1 



where 



r,q,m,n=—M/2 
M/2-1 



r=-M/2 
M/2-1 

m^n=-M/2 J 
( M/2-1 M/2-1 

y r=-M/2 q=-M/2 
M/2-1 1 

+ Y Fm,mFlnW{X'Z^, ^"WCKL ^''') \ + 
m,n=-M/2 J 

f M/2-1 _ ^ ^ . . 1 

\^m,n=-M/2 J 



M2 

1 

M2 



J\//2-l 
n=-A//2 



(0^ 



M/2-1 



ej^r=-M/2 

Using the results of fl5.2ip . we can express the first term of fl5.12p as follows 



1^ 



(01 
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-E 



^(0 



00) 



dt 



M\e 



N/2-1 
n=-N/2 



(01 



M2 



M/2-1 
n=-M/2 



M 



|^(0|2^(0 
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Using some of the arguments we used to reduce the first term of fl5.12p . the second 
term of fl5.12p can be simphfied as follows 



2\9 



(01 



M 



E|y"Re|^rf^''(t)rf«*(t)| dt 



Qm\ 



M 



-Re < E 



^(0 

'em 



M 



Re<^ E 



M J 

7V/2-1 



Ei / S'\t)d'^^>{t)dt 



A//2-1 



M2 



M2 



Re 



n=-N/2 
N/2-1 

n=-N/2 



m=~M/2 
2 - 



M2 



N/2-1 
n=-N/2 



(5.25) 
(5.26) 
(5.27) 

(5.28) 

(5.29) 



Equation fl5.27p follows from fl5.15p . (15.170 . and the facts that E{6„6j^} = 6n-m 
and E{c„cJ^} = Sn-m- Finally, the third expectation of (I5.12p can be reduced as 
follows 

|^(0|2 



M 



(5.30) 



Thus, combining the results of (I5.24p . (15.290 . and (I5.30p . we can express (15.121) as 



M 



-E 



e^'^S'\t)-e^'^d^'\t) 



M2 



M/2-1 
n=-M/2 



(0^ 



(5.31) 



Now by combining the results of (15.111) and (I5.3ip . -gE |||ey(t)||^} can be upper 
bounded as follows 

L 



>{iM*)in<Eri^-^ 

1=1 



M/2-1 

J2 F„,„iy(A^:;„,/3 

n=-M/2 



|0"f cD«. (5.32) 
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By applying the results of Appendix O and Appendix [Hi we find 

^.E{||e,(t)||^} 



M 



-M 4 '''^^ 



=1 

X I smc I I ^^-r^ I i + z I - z I — I 











{ 1 - a(0 J 


' 1 - a(') / 



+ + . ,5.33) 

where B < n is chosen to make sinc(a;) > cos(i?x) for |x| < ^. In Section [5.2.1^ we 
examine how tight this bound through simulation. 

5.2.1 Simulation of Estimation Error Using Doppler Approx. 

To examine the tightness of the bound in (15.331) . we ran a simulation of a one 
path channel similar to the simulation described in Section 14.3.11 The single path 
channel again had parameters given by a'^' = —0.001, and 7*^' = 0. A training 
signal modulated of 300 random symbols was sent through the channel. This train- 
ing signal had a fractional bandwidth jj; ~ ^- Parameter estimates were given by 
^(1) _ ^(1) _j_ ^ ^j^^ g(i) _ ^(1) -Jf-^, Thus, we simulated the worst case scenario where 
the parameter estimates were half a grid spacing away from the true parameter val- 
ues. The estimate of 6^^^ was made based on (15. Sp . Next a data signal modulated 
from a series of 300 random data symbols was sent through the channel. The error 
metric -^E {||ey(t)||^} was calculated for different many choices of grid spacing. This 
simulation was repeated 10 times and the results were averaged. We set B = 1.8138 
and plugged in the simulations parameters into (15.331) to calculate the error upper 
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bound as 

< 1.004- 0.2492 cos2(^:5569^) 

X (sine (149.925lAa + 0.0865) + sine (149.9251Aa- 0.0865))^ (5.34) 

Figure 15.21 plots this upper bound while Figure 15.11 plots the average simulated 
error. Figure 15.31 shows the two together, and it can be seen that the simulated error 
is less than the upper bound we derived. These results are also shown in Figure 15.41 
and Figure 15. 5[ Figure 15.41 plots the simulated error and the derived bound for fixed 
values of Aq, while Figure 15.51 plots the simulated error and the derived bound for 
fixed values of A^. It all cases it can be seen that the bound derived in (15.33^ does 
indeed hold and appears to be relatively tight. 

5.3 Doppler-shift Model vs. Scale-shift Model 

A primary motivation of this thesis is to analyze the cost of using the Doppler 
approximation. In Chapter H] we derived an expression fl4.67p which upper bounds the 
average estimation error when a model composed of the effects scale-shifts is used. 
Then in Chapter [S] we derived a similar expression (15.331) which upper bounds the 
average estimation error when a model composed of the effects of Doppler-shifts is 
used. Through simulation, we have shown that both these upper bounds are relatively 
tight. From these simulations, we can also see that there is an additional cost occurred 
when using the Doppler approximation. For instance, by comparing Figure with 
Figure [57il we see that when A^ = A^ = 0, the simulated error is when the scale- 
shift model is used, but when a Doppler-shift model is used, the simulated error is 
roughly 0.024. 
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A„ X 103 



Figure 5.1: Average simulated error using a Doppler-shift model 
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Figure 5.2: Derived upper bound for a Doppler-shift model 
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Figure 5.3: Derived upper bound vs. simulated error for a Doppler-shift model 
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Figure 5.4: Simulated error and derived upper bound vs. for fixed values of A^ 
using time-sliift/Doppler-shift model. Simulated error is plotted with solid lines while 
the derived bound is plotted with dotted lines. 
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Figure 5.5: Simulated error and derived upper bound vs. for fixed values of ^ 
using time-shift/Doppler-shift model. Simulated error is plotted with solid lines while 
the derived bound is plotted with dotted lines. 
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We now analyze the effects of the fractional bandwidth and rate of change of 
the channel paths {a'*'};^^ on the error resulting from the Doppler approximation. 
To do this, we consider the case where L = 1 and M ^ 1 to approximate the upper 
bound of (I4.67P as 



wideband 



(0|2 



1-2011 


-^|^"f cos^j' 


M ' ' 





,(0 



2T 

B \ MA„ 



1 - aC) J 477 
B \ MA„ 



and the upper bound of (15.331) as 

1.2011 



Doppler 
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^(0|2_ ±1^(012 pQg2 



1 - a(" / An 



fBA^ 



(5.35) 
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sine ( nf/rfl H ^ 

sine ^vr/cT^^l + 



\ 2T 

1 \MAa B Ma^'^ 



Ait 
MA„ 



B Ma^'^ 



n 2 



. (5.36) 



ad) J 47r l-a« 2n 

By considering the difference between £^widoband and £^Doppior for practical (relatively 
small) values of £^widoband, we can analyze the cost of the Doppler approximation. 
Notice that the scale-error contribution to £^wideband is proportional to the magnitude 
of 



•'wideband 



( 7r/eTfl + i \ +- 



B \ MAn 



a'*' / An 



(5.37) 



In other words, when |0wideband| is very small, the scale-error contribution to £^wideband 
is also very small. Similarly, the scale-error contribution to £^Dopper is proportional to 
the magnitude of 

MAa B Ma'-'^ 



•'Doppler 



vr/cT(l + ^ 
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An 



l-a« 2n 



(5.3^ 



Now we consider a fixed value of </)„idcband and solve for (^ooppier in terms of </>wideband- 



Solving for in (15.371) and plugging the result into (I5.38p . we can write 

, _ , _ B [ Ma^'^ (/>widcband 

vDopplcr *Pwidcband i I n / \ 

From (I5.39p . we notice that as fcT — >■ oo, we get 



(5.39) 



B Ma . 

0Doppler — ^wideband ^ • \pA\}) 

i — a ZTT 

This shows that there two conditions in which the Doppler approximation is accurate. 
One is that the signal should have a small fractional bandwidth (i.e. large fcT). The 
other is that a < 0widcbandj^ = 4>wideha.nd^ (^mt)1(i/t) where MT can be recognized 
as the signal duration and 1/T can be recognized as the signal bandwidth. 

To further illustrate the effects of fcT on the error resulting from the Doppler 
approximation, we plot the upper bounds given in (15.351) and (15.361) over many values 
of fcT in Figure EH For this plot, we set L = 1, a^'^ = .0008, B = 1.8128, M = 100, 

= 0, and chose so that i^widcband was constant for all fcT. As can be seen, the 
approximation becomes more accurate as fcT increases. 
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Figure 5.6: Upper bound of (14.671) and upper bound of fl5.33p vs. fcT 
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CHAPTER 6 



Experimental Results and Conclusion 

6.1 Simulated Channel 

For further analysis, we simulate the estimation strategy described in Chapter [3] 
over a three path channel. In this simulation, the training signal we use has a carrier 
frequency of 10 kHz, a baud rate of 10 kHz, and a sampling frequency of 30 kHz. 
The transmitted training signal is modulated from 200 random QPSK symbols. In 
order to estimate the channel parameters, an under-determined Matrix composed of 
the effects of time-shifts and scale-shifts on the training signal is constructed. We 
simulate the worst case scenario where the true channel parameters are half a grid 
spacing away from the parameter estimates (i.e. a*'' = a*'' + ^ and 7*'' = 7''' + ^ 
for I = 1,..,L). Figure WTl and Figure \Q2\ show the true channel parameters and 
parameter estimates found using Orthogonal Matching Pursuit (OMP) [5]. 

To analyze the estimation strategy further, we define two error metrics. The first 

II " 1 1 2 

z— A(? 

we call the training error which we define as " a " , where z is a vector of noisy 

1^1 

measurements of the received training signal, A is a matrix composed of the effects 
of time-shifts and scale-shifts of the training signal, and is the sparse estimate 
of the channel gains. The second we call the the data error which we define as 

54 



S 0.5 - 
1 



O Actual Channel Parameters 
Strongest Estimated Taps 



11.5 12 12.5 13 13.5 
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Figure 6.1: Plot of true channel parameters (red circles) and the grid restricted 
parameter estimates (black diamonds). The bottom axis ranges from 7=10ms to 
7=15 ms. The left axis ranges from a=-2.5 ms/s to a=2.5 ms/s. Figure ESI shows 
another view of this channel estimate. 
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Channel Estimate 
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Figure 6.2: Another view of the channel estimate shown in Fig. 16.1! 



; where y is the received data signal (noiseless), and the columns of B are 
composed of the effects of time-shifts and scale-shifts on transmitted data signal. 
In this simulation, the channel seen as the data signal is sent is the same as the 
channel seen as the training signal is sent. However, we know that algorithms, such 
as least squares, can minimize the training error with inaccurate (noise corrupted) 
estimates [3] . For this reason, we are most interested in the data error for our analysis. 

Figure 16.3! shows both the training error and the data error vs the number of 
iterations using OMP for the noiseless case. We see that both the training error 
and the data error decrease as the number of OMP iterations increases. However, 
when the received training signal z is noisy, we find this is not the case. Figure 16. 4[ 
Figure [63| and Figure [6l6] plot the results for different signal-to- noise ratios (SNRs). 
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As can be seen, the data error begins to increase as the number of OMP iterations 
increases. The lower the SNR, the smaller the number of OMP iterations which 
minimizes the data error. This illustrates how a sparse solution (few OMP iterations) 
can give a better solution which is less corrupted by noise. 



Training Error 
Data Error 




10 20 30 40 50 60 70 
OMP iterations 



80 90 100 



Figure 6.3: Estimation error vs OMP iterations (noiseless case). 



In addition to OMP, we also test Basis Pursuit (BP) [7] which finds the solution 



to 



minA|0|i + 1/2 
e 



z-AO 



For this algorithm, we measure the data error and the training error over several 
values A. Figure 16.71 Figure 16.81 Figure 16.91 Figure 16.101 plot the results for four 
different SNRs. From these figures we see that the lower the SNR, the higher the 
value of A which minimizes the data error. 
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Figure 6.4: Estimation error vs OMP iterations (SNR = 1 dB). 




-2 
-4 
-6 
dB -8^ 
-10 
-12 
-14 
-16 




Training Error 
Data Error 



10 20 30 40 50 60 70 80 90 100 
OMP iterations 



Figure 6.5: Estimation error vs OMP iterations (SNR = 3 dB). 
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Figure 6.6: Estimation error vs OMP iterations (SNR = 10 dB). 




Figure 6.7: Estimation error occurred using Basis Pursuit (noiseless case). 
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Figure 6.8: Estimation error occurred using Basis Pursuit (SNR = 1 dB). 




Figure 6.9: Estimation error occurred using Basis Pursuit (SNR = 3 dB). 
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Figure 6.10: Estimation error occurred using Basis Pursuit (SNR = 10 dB). 

Finally we test Fast Bayesian Matching Pursuit (FBMP) [8]. This algorithm 
has a user defined input p which corresponds to the probability of a non-zero tap. 
Figure I6.1H Figure 16.121 Figure I6.13[ and Figure 16.141 plot training error and data 
error vs p for different SNRs. We find that FBMP tends to outperform both OMP 
and BP for this application. 

6.2 Unknown Channel 

To study the estimation strategy further, we implemented it on experimental data 
collected from the Woods Hole Oceanographic Institution. Specifically, we used a 
transmitted passband signal and a received passband signal to estimate an unknown 
channel. The transmitted signal had a carrier frequency of 13 kHz, a bandwidth 
of 10 kHz, and a time duration of about 60 seconds. Both the transmitted signal 
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Figure 6.11: Estimation error occurred using FBMP (noiseless case). 
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Figure 6.12: Estimation error occurred using FBMP (SNR = IdB). 
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Figure 6.13: Estimation error occurred using FBMP (SNR = 3 dB). 
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Figure 6.14: Estimation error occurred using FBMP (SNR = 10 dB). 
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and the received signal were sampled a frequency of 39062.5 Hz (roughly 4 times the 
signal bandwidth). For the first channel estimate, only the effects of time-shifts were 
considered as done in Section 14. 2[ Estimates were repeatedly made using OMP over 
several time periods. Fig. 16.151 and Fig. 16.161 plot the resulting estimate . As can be 
seen, the channel is multi path and varies with time. 
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Figure 6.15: Channel Estimate using OMP and a LTI model. Bottom axis represents 
time in seconds and the left axis represents the channel delay in milliseconds 

Since the true channel parameters are of course unknown, we can not use the 
same error metrics we used in previous sections. For this reason, we us the following 
metric, similar to one used in [2]. We denote this metric as the prediction error and 
define it as follows 
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Figure 6.16: Estimate shown in Fig. l6.15l on log scale 
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el[i + l] = \z[t + 1] - z[i + 



(6.1) 



where z[2 + 1 \i] is the prediction of the received signal at time i + 1 based on the 
signal received up to time i. 

Figure 16.171 plots this prediction error for various grid spacings. For this simula- 
tion, a LTI invariant model was constructed to estimate 201*^ sample ^[201] based on 
the previous 200 samples z[l : 200]. This process was repeated until estimates were 
made for z[201 : 300]. We normalized the prediction error so that Figure 16.171 plots 
||z[201 : 300] - z[201 : 300]||V lk[201 : 300] ||^ vs ^. As can be seen, the prediction 
error increases with increasing grid spacing. 




Figure 6.17: Normalized prediction error vs. ^ 
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6.3 Conclusion 



In this thesis, we have analyzed the use of spare reconstruction algorithms for 
the application of channel estimation. In particular, we examined the effects of the 
choice of grid spacing on the average estimation error. We derived expressions which 
upper bound the average estimation error in terms of the choice of grid spacing for a 
a linear time invariant model, a time varying narrowband model, and a time varying 
wideband model. Additionally, we explored the cost of using a narrow band model 
to approximate a wide band model. We then quantified the cost of this approxi- 
mation in terms of the fractional bandwidth of the signal and the rate of change of 
channel. Finally, we compared through simulation the use of three different sparse 
reconstruction algorithms (OMP, BP, and FBMP) for this apphcation. 
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APPENDIX A 



Expected Value Evaluation 



In this appendix, we simplify -^.E \ 9<-'^ d'-'^ (t) - O'^'W^'^t) 



as follows 
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(A.2) 



(A.3) 



Equation ([QD can be reduced further by noticing that the sum of the second term 
gc/«(t)rfW*(t) and third term ^S^>{t)d^'^{t) is equal to 2Re ||j}c;(')(t)rf«*(t)}. 
This is because for any G,F eC, G*F + GF* = 2Re {GF*}. We thus rewrite ([Ql) 
as 



1e 

M 



1^ 



(01 



M 



-E 



^(0 



^(0 



d^'\t) ' - 2Re <! ^S'\t)d^'^*{t) \ + \d^'\t)f dt \ (A.4) 
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Equation flA.4p can be split into 3 different terms as follows 

2 



M 



(01 



M 



-E 



2\9 



(01 



M 



-E< / Re 



eh 



dt 



em 



S'\t)d^^>{t) ) dt 



(01 



M 



-E 



\d''\t)\^ dt 



(A.5) 
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APPENDIX B 
Upper bound on a 



In this appendix, we derive an expression defined in (14.260 . in terms of A^. 
To do tliis, we first define 5''' = ^^{7''' ~ 7*'^}- We assume that the — 7*''| < j, 
so that l^^'^l < 0.25. Now we use the change of variables m = q — r to write 

"^'^ = ]V2 5Z {N-\m\) R\mT + 6^'^T). (B.l) 

m=-N+l,ni=iO 

Defining has the roll off factor of the raise-cosine pulse R, we can write 



R{xT) 



cosi^irajix) sin(7rx) 
1 — {2aRxy TTx 



(B.2) 



By noticing that 



cos{TTafix) 



1 - i2anxy 



< 1 for all aji, x 



(B.3) 



and that | sin(7rx)| < 1, we can write 



R^{xT) < 



(B.4) 
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Putting fIB.ll) together with flB.4p . we get 

N-l 



a 



(0 < 



1 



N — \m\ 



1 ^-1 



N — m 



N-l 



+ 



—y 

2/V2 



N — m 



< 



m=l ^ ' m=l ^ ' 

2 iV-m 

^2^2 



(m - 0.25)2 



2N -0.5 



E 



iV-0.25 
(m- 0.25)2 



Af-1 



E 



m - 0.25 
(m- 0.25)2 



E 



E — 



(B.5) 
(B.6) 
(B.7) 
(B.8) 
(B.9) 



vr2Ar2 0.25)2 7r2iV2 ^ m-0.25' 

r?i=l ^ ' m=l 

where f lB.7p is due to < 0.25. We can upper bound the first term in (IB. 90 as 



follows: 



N-l 



E 



16 



(m - 0.25)2 



7V-1 

+ E 

m=2 

N-l 



>- 0.25)2 



< 



16 ^ 1 
"q" + 2^ (m- 1)2 



9 

16 



m=2 
N-2 



9 ^ m2 

m=l 
16 7T^ 



using the fact that J2m=i ~ \ ^^"^ noting that all the terms in the series Yl'^=i 



(B.IO) 
(B.ll) 

(B.12) 
(B.13) 

oo 1 



are positive. To upper bound the second term in flB.9p . we lower bound its negative: 



N-l ^ N-l ^ 

y — - — > Y- 

^ m - 0.25 ^ rr. 



m=l 



m 

m=l 



(B.14) 
(B.15) 
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Plugging flRTSjl and flRlSjl into flROD . we get 

2iV- 0.5,16 7r\ 21n(A^) , 
694 

< (B.17) 

Note that, when is large, 2A^-0.5 ^ 2N and (f + 7rV6)iV-ln(A^) ^ (f + 7rV6)A^. 
Thus, although the bound in flB.161) is a tiny bit more accurate than the one in flB.17p . 
the latter is much more intuitive. 
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APPENDIX C 



Inner Product Evaluation 



To simplify the inner product we first express it using ( 13.29^ . (13.2ip . 

and ([322]). 

s'"(t)#*(t)dt (C.l) 

^i27r/c({l-a(0)7(0-(l_a(0)7(0+{a(0_a(0)t) 



M/2-1 

X Vl - a« ^^^((^ " - rT - (1 - a^'^)Y'^) 

r=-M/2 
M/2-1 

X VT^r^ J2 hlp{{l - a^'')t - qT - {I - a^'^)f'') dt (C.2) 

q=-M/2 

M/2-1 M/2-1 

= ^/r3^^/r3ai^e^-2-/c((i-(0)^(0_(i_.(0),(0) ^ ^ 

r=-M/2 q=-M/2 

X y e^^2-/=('^'"-'^*")V((l - ci^'')t - rT - (1 - a«)7*'0 

X - a('))t - gT - (1 - a('')7«) dt (C.3) 

To evaluate ( 1C.3I) . it is helpful to approximate the term e-'^'^'^'^*^'^*''' with a time 
invariant constant in order to evaluate the integral. To do this, we note that the 
term — d*''))f: — gT — (1 — a''')7''') has its maximum value when t = ^^rfpy + 7*'\ 
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but diminishes to as |(1 — a*'')t — gT — (1 — a*'')7*''| increases. Likewise, the term 

— a*'')t — rT — (1 — a*'')7*'') has its maximum value when t = + 7*'' 

diminishes to as |(1 — a*'')t — rT — (1 — a*'')7*''| increases. We define F^'^, as the 
value of e^'^'^fcia'-'-^-'i'-'^}t ^vhen t is half way between the peak values ti = x^fpy + 7*'' 
and t,^j^ + Y'K 

Fi^l A (C.4) 
The term eJ'2^/-("^'' in ([03]) can be well approximated by F^'^l over the periods of 



t where both - a(")t - gT - (1 - a«)7«) and - a('')t - rT - (1 - a('')7(')) are 
not negligible. This in part is because we make the mild assumptions that the value 
of a^'^ to be relatively small (i.e. a*'' < 0.01) and that d*'' ~ a*''. It is also because 
we have assumed that the signal has a high fractional bandwidth (i.e. ~ !)• By 
using F^'^, to approximate e-?^'^-/^'=("*''~"*'')*, we can express (1C.3P as follows 

(s<')(t),S<"(t)) 

M/2-1 M/2-1 

= ^/r3^^/r3ai^e^-2-/c((i-a(0)^(0-(i-,(0),(0) ^ ^ 

r=-M/2 q=-M/2 

X H'!r j - «^")^ - rT - (1 - a('07«) 

X p((l - a^'^t - gT - (1 - a«)7«) dt (C.5) 



Now if we make the substitution v = {1 — d''')t — gT — (1 — a''')7*'', (IC.5I) becomes 

/- A//2-1 Af/2-1 

Vi^a^^j2^/,({l~a(0)7(0_(i_a(0)T,(0) 6*T'" 

^ r=~M/2 q=-M/2 

X / P{v)p{--—v + — - rT + (1 - a''))(7<'' - 7*'')) 

A//2-1 A//2-1 

= ^/3(0e^-2-/c((l-<")7«-(l-a(0),(0) ^ ^ ^^^*^(o 

r=-M/2 q=-M/2 

X y p(t;)p(/3'"t; + /3<"gT-rT+ (1 -a'")(7'" -7<'0)dw, (C.6) 
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where = ■ To simplify flC.6p further, we use the Wide Band Ambiguity 
function defined in [13] 



W{X, /3)^^P p{t)pm + A)) dt. (C.7) 



Using this function, we can express flC.6l) as 

M/2-l M/2-1 
r=-M/2q=-M/2 

rji 

xW{qT - ^ + (1 - a'")(7('> - t'"),/?*")- (Ci 
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APPENDIX D 
Upper bound on Cj 



In this appendix, we derive an expression ti)''' (defined in (I4.54p ) in terms of M . 
To do tliis, we first note that the wideband ambiguity function (defined in (1C.7I) ) can 
be simplified as follows 



W{\(3) ^ ^ (3 p{t)p{(3{t + \))dt (D.l) 



j j P{h)e^^^f^'dh j P{h)e^^''f'^^'^''^df2dt (D.2) 
-^11 P{h)e^'-f^'dh j P{f/(3)e='-f^'+'^^dfdt (D.3) 
j P{h)P{f/(3)e^^^f^ j e^^-^^^+^^'dtdfdh (D.4) 



s{h+f) 



j P{-f)P{f/P)e^'-f^df. (D.5) 



We assume that p{t) is the SRRC pulse with rolloff = so that 
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in which case 



min(l,3) 

^ r 2T 



W(A,/3) = -= / e^'-f^df (D.7) 



^ "^^"^^'^)sinc(^^^^^A) (D.8) 



^ T ' T 

ii(i 

T 



T min(l,/5)sin(7rHH2i^A) 



rp2 



Using fID.lip . the term ci)''' can be upper bounded as follows 

M/2-1 A//2-1 



■,(0 A 

r=-M/2 e^r=-A//2 



(D.9) 



sin(7r ^T^A) (D.IO) 



W(XJ3Y < ^ (D.U) 



M/2-1 A//2-1 

-A-//2 e^r=-A//2 ^ ^e,r P 

We assume that a*'^ < ^jf and a*'' < so that for large M, we can make the 
approximation = \zt(i) ~ 1 + — a'''. Using this approximation, we can express 
A^''^ (defined in flOOD ) as follows 

A« = rT - eT + (a'" - a'^^eT - (1 - a«)T(5''> (D.14) 
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where 5'-'' = ^{^'^ ~7'''} ^ j- Using the change of variables m = e — r and combining 
flDl3|) and fiDlil) . it follows that 

M/2-1 M/2-1 

M-l M/2-1 

= t y y I 

-1 Af/2+m 

+ ^-- E E 7 ^ ^ (D-16) 



A/-1 

Z 

< 

M%2/3(o + (^(0 _ a(n)e + (1 - a(o)5(o)' 



A/-1 Af/2^1 



=1 e=-A//2 

Notice that, since |e| < ^, 5''' < i and a*'' < then 



|(a«-a'")e + (l-a«)5'''| < | (a<" - a''')e| + | (1 - a''^^''' | (D.18) 

= |a'"-a«| |e| + |l-a<"| (D.19) 

< i«"-«">if + a + i'i 

We assume that |a^'' — a*'^| < — With this assumption in mind, it imme- 
diately follows that |(d*'' — a*'')e + (1 — d''')(5'''| < ^. In this case 

Af-l 

cl;"' < = — ^ (D.22) 

~ M7r2/?(0 (m - 0.5)2 

A/-1 \ 



M7r2/3(0 \ ±^Am- 0.5)2 



2 

< 



m=2 
A/-2 



, 4 + E A 1 (D-24) 

2(4 + 7rV6) 

< ^ / ^ p.25 

M7r2/5w 
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^ 1 TT^ 

using the fact that — - < —. Finally, since 

1 1 - a« 



we have 



6 

m=l 



^(0 1 - a(') 



< l + |a(')-a(')| 

< 1 + -(1--) 

1 2M + 1 

< 1 + = —. 

2M 2M ' 



2(4 + 7r76) 2M + 1 



Mtt^ 2M 

Now if we say that M > 10, we know that < §, so that 

.(0 < 2(4 + 7rV6) 21 ^ 1.2011 

Mtt^ 20 - M ■ 
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APPENDIX E 



Lower bound on 



M/2-1 



n=-M/2 



The primary goal of Section 14.31 is to find an insightful expression which relates 
an upper bound of the error metric -pE 1 1 \ey{t) 1 with the the grid spacing (A^ and 
A^). With this goal in mind, it is necessary to find a lower bound on the expression 

M/2-1 ^ 
n=-M/2 



0^ 



that I a''' I < ^ and |a<'> 



in terms of the grid spacing. For our analysis, we assume 

— 2l?- assume that the bounds on the quantities 

— aj^d ~7'''| are slightly tighter than the bounds assumed in Appendix [D| 
so that 



2AI\^ 2M)2fcT 



(E.l) 
(E.2) 



Note the presence of the fcT term in (IE. 11) . and recall that fc>^ (assuming that 
the carrier frequency is centered in the signal passband) which implies that 2/cT > 1. 

We start by recalling the definitions of and W^(A,5^'„, /^''') given in (]C.4p and 
(1C.7I1 . respectively. We now notice that 



n=-M/2 



M/2-1 
n=-M/2 



(E.3) 
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where 



X^' = vr/,(a<')-a''))(^ + ^). 
Using the Euler expansion K^^ = cos(x*^') + j sin(x^'), we get 

2 



(E.4) 
(E.5) 



M/2-1 
n=-M/2 



{^)^ 



M/2-1 

E W^(^^!n,/^<'')[cos(x^:0+jsin(x^:0] 

n=-M/2 

(M/2-1 ^ ^ 

E W^(^n!n,/3<'')C0S(XJ:'^ 
n=-M/2 



(E.6) 



(E.7) 



where the inequahty follows because iy(A^^„, G M. The latter bound is actually 
quite tight because the quantity W^(A*^'„, is symmetric around n = 0, the quantity 
sin(x^^) is antisymmetric around n = 0, and the summation is almost symmetric 
around n = 0. 

Next we note that, under our assumptions, we have 



I ^n,n I 



Y^nT{l - a« - 1 + a'')) + (1 - a''))(7« - 7 



(0^ 



\n\T _|_ T ^ T 
2M-2fcT 4 — 2 



1 



2M ^ 2/cT 



^ + (i + 5I7)-t(i + 



2M' 



lAn I — "Jc 2MV 2M>2fcT 1-^ 2M 



7r|n| 



(E.8) 
(E.9) 
(E.IO) 
(E.ll) 
(E.12) 



Due to the range of |Xn'l? term cos(x^') is clearly non-negative over the range 



of n in (Ell). We now show that W {X'i]^, p^'^) 
non- negative over those n. First, notice that 



miii(l,^W) 



mm 



> 



iin(l,4W) 



1-- 
I 

1+1 



smc 



iin(l,/3W) 
T 



M-1 
M+1 



K]n) is also 



(E.13) 



because - — ^ < < ¥ . Then, notice that the argument of the sine obeys 

Si^^A^^) < ^ < i (E.14) 



Because both cos{Xn) and Vr(A^'„, are non-negative over the range of n in ( 1E.7P 



we can lower-bound (1E.7P by plugging in lower bounds for W^Xl^^^, in the sum- 
mation. Next, we derive this bound. 

Because sinc(a;) decreases as x increases for |a;| < | 



W^(A^:!„,/3''')>y^sinc(4t). (E.15) 



Applying this inequality. 



M/2-1 
n=-M/2 

W^AY. sine (f^) cos (^vr/,r(a« - a<0),) . (e.16) 

\ n=-M/2 / 

We now observe that for particular positive values of i? < vr, sinc(x) > cos(i?x) for 
< I (see Fig. lE.ip . Since < |, we can rewrite (1E.16P as 



M/2-1 
^n=-M/2 

, / M/2-1 
M-1/ V- /baW 



^ mTI E COS (^) COS (.r/.T(a<o - a<> (^(0 + 

/ A//2-1 

T E cos(^n(a«-a(')) + f(l-a«)(fO_^(0)) 



M + 
M-1 



M + 



, n=-M/2 



X COS ( 7r/.T(a'0 - a<'))n ( + ) ) ) (E.17) 



M - 1 



, 2 

A//2-1 



M 



- E cos (^C'^'n + cos (^E«nH , (E.18) 



^n=-M/2 
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Figure E.l: Plot of sinc(x) and cos(-Bx) when B = 1.8138 



where 



D"^ ^ ^(l-a''')(f"-7('0 (E.20) 
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Equation flE.lSp can be recognized as a Riemann-sum approximation of the following 
continuous integral 

M-i ( r^'^ 



M +1 



-Ml 



COS (c^H + cos dt 



M-l [^''^ ( ( \ ( 

- sin [C'^'^t^ sin (^*")) cos (^E<"t) dt^ (E.22) 
= ^ y ^ cos [C^'^t^ cos cos dt\^ (E.23) 

Equation flE.22p follows from the trig identity cos(x+?/) = cos(x) cos(?/) — sin(a;) sin(?/) 
and flE:23|) follows from the fact that /^(^^ sin (^C^'^t) cos (E^^'^t^ dt = 0. We can now 
use the trig identity cos(x) cos(?/) = | cos(x + y) + | cos(x — y) to rewrite ( ]E.23p as 
follows 

M-l f r''/^ 

'-A//2 



M+ 1 

M-l 



M + 

M-l 
M+ 1 

M-l 



/" cos cos cos (^E^'H^ dt 

^ cos2(D«) ^ i COS + + i COS - j 



COs2(Z) 



(0^ 



£■(0+0(0 ' £(0-c(o I (E.24) 



^ ^ - cos2(D«) sine + C«)f ) + f sine - C('')f )) (E.25) 

■ ^ cos2(D«) (^sinc + ) + sine (^(E"' - C«)f ))|e.26) 



84 



Finally, by combining flE.7p . flE.lSp . and flE.26p and plugging back in defined expres- 
sions for C^'\D^'\ and E^'^ (equations flKlOl) . ^K20\ . and ^K21^ ). we find 



M/2-1 
n=-M/2 

M-1 M2 



> 



M+1 4 
X I sine 



cos2(f(l-aW)(7W-7«)) 



+ 



+ sine 



1 - a(') 1 - y 27r 



1 - a(') 



1 - a(') 



277 



(E.27) 



Since cos(x) and sinc(x) are symmetric around x = and we have assumed |a'''| < 



q(i) _ q(0| < g^j^^ 1^(0 _ ^(o| < fii^ gg^j^ j^Q^ write 



M/2-1 
n=-M/2 

> ^.^cos^(^(l + ^)) 

7r/,T , 7rf,T + B\ MA. 

47r 



X smc 



1 - 



+ 



2M 



1 - a('' 



^ . , , 7r/,T nf,T-B \ MA, 

+ smc ^ — 

1-217 / 4vr 



(E.28) 
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APPENDIX F 



Inner Product Evaluation (Doppler Case) 

To simplify the inner product {s^'\t), s^'\t)), we first express it using ( 13.29^ . (13.2ip . 
and (O). 

= J gi2^/c((l-a(0)7{0-(l-a('))T,(0 + (a(0_a(0)t)^/Y3^ 

M/2-1 M/2-1 

X J2 brp{{l - a^'^)t - rT - {1 - a^'^h^'^) ^ b^pit - qT ~ ^^'^) dt 

r=-M/2 q=-M/2 

M/2-1 M/2-1 

r=-M/2q=-M/2 

^^((1 - a^'^)t - rT - (1 - a'")7<")p(t - - 7'") tit (F.l) 

To further simphfy (IF.ip . it is helpful to approximate the term e-^^'^-^'^*^"'''""*'^)* with a 
time invariant constant in order to evaluate the integral. To do this, well note that the 
term ^(t — gT — 7*'') has its maximum value when t = gT + 7*'', but diminishes to as 
|t — gT — 7*''| increases. Likewise, the termp((l— a*'')t— rT— (1— a*'')7''') has its maxi- 
mum value when t = y3^+7*'' and diminishes to as |(1 — a''')t — rT — (1 — a''')7'''| 
increases. We define F^'^l as the value of e-'^'^-'^'^*^''*''""*'')* when t is half way between 
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the peak values ti = qT + 7''' and t2 = j^^o + 

The term eJ'2'^^-(^*''-'^''')* in (D) can be well approximated by F^'J. over the periods of 
t where both p{t — qT — 7*'') and p((l — a^'^)t — rT — (1 — a''')7*'') are not negligible. 
This in part is because we make the mild assumptions that the value of a^'^ to be 
relatively small (i.e. a^'^ < 0.01) and that d^'^ ~ a^'K It is also because we have 
assumed that the signal has a high fractional bandwidth (i.e. ^ 1). By using F^'-J. 
to approximate e^'^'^-^"^"'^'^'"'^'^^^ , we can express (IF. 10 as follows 

A//2-1 M/2-1 

r=-M/2 g=-M/2 

X Fg^r j p((l - - rT - (1 - a<'')7*")p(^ - qT - 7«) dt. (F.3) 

Now if we make the substitution v = t — qT — 7^'', (IF. 30 becomes 

M/2-1 A//2-1 

^/^3^^,2./.({l-a(0)^(0_(i_„(0),(0) ^ ^ ^^^.^^^^ 

r=-M/2 g=-M/2 



X / p(t;)p((l - a^''>)v + (1 - a«)gT - rT + (1 - af")(7« - 7^)) dv 

M/2-1 M/2-1 

^(,)gi2^/c((l-a(0)^(0_(i_a(0)^(0) ^ ^ ^ 

r=-M/2 q=-M/2 

X y p(t;)p(;3<S + ;9'"gT-rT + /3<'>(7<" -7<'>))rft;, (F.4) 
where = 1 — a^''. To simplify (1F.4P further, we use the wide band ambiguity 



function defined in (10.70 . so that we can express (1F.4P as 



M/2-1 M/2-1 

r=-M/2 q=-M/2 

xPF(gT-^j + f')-7«,/3<"). (F.5) 
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APPENDIX G 



Lower bound on 



M/2-l 



i=-M/2 



The primary goal of Section 15.21 is to find an insightful expression which relates 
an upper bound of the error metric -pE 1 1 \ey{t) 1 with the the grid spacing (A^ and 
A^). With this goal in mind, it is necessary to find a lower bound on the expression 

M/2-1 ^ 
n=-M/2 



that I a''' I < ^ and |a<'> 



in terms of the grid spacing. For our analysis, we assume 

— 2l?- assume that the bounds on the quantities 

— aj^d ~7'''| are slightly tighter than the bounds assumed in Appendix [D| 
so that 



|;y(0_^(0| < 



2M\^ 2M)2fcT 

T 2M-2 
4 2M-1' 



(G.l) 
(G.2) 



Note the presence of the fcT term in (IG.ll) . and recall that fc>^ (assuming that 
the carrier frequency is centered in the signal passband) which implies that 2/cT > 1. 

We start by recalling the definitions of and Vr(A^'„, given in (lF.2p and 
(IC.7I1 . respectively. We now notice that 



A//2-1 

n=-M/2 



M/2-1 
n=-M/2 



(G.3) 



where 



-(0 



^(0 A ^3Xn 

7r/,(5(')-a''))(nT + ^). 



yd) A 



(G.4) 
(G.5) 



Using the Euler expansion K^^ = cos(Xn) + J sin(x^'), we get 

2 



M/2-1 
n=-M/2 



M/2-1 

E W^(^^!n'/^^')K(X^:')+Jsin(x«)] 

n=~M/2 

(\ 2 
M/2-1 
n=-M/2 



(G.6) 



(G.7) 



where the inequahty follows because W^(A*^'„, G M. The latter bound is actually 
quite tight because the quantity Vr(A^'„, is symmetric around n = 0, the quantity 
sin(x*^') is antisymmetric around n = 0, and the summation is almost symmetric 
around n = 0. 

Next we note that, under our assumptions, we have 



m = lTzi(i)^r(l-a<')-l) + (fO_^ 



2Af-l 



< ? 



(G.8) 
(G.9) 
(G.IO) 



< ^/c • 217(1 -2i?)2iT- Mm + IZV) < |H (G.ll) 



< 



4 ' 



(G.12) 



where flCTTOl) and flGa2l) follow from the fact that \n\ < f. Due to the range of 
Ix^'l, the term cos(x[J') is clearly non-negative over the range of n in (]G.7p . We now 
show that Vr(A[J|„, = sinc C^™*-;^^' ''' ^n|n) ^^^^ non-negative over those 

n. First, notice that 

^ ' (G.13) 



n(l,/3W) 



mm 



J 
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> 



M 



because 1 — jj < < 1 + Then, notice that the argument of the sine obeys 

"^^^¥^>^Tn < # < i (G.14) 



Because both cos(x^O and (A[J^^, are non-negative over the range of n in flG.7p 



we can lower-bound flG.7p by plugging in a lower bound for Ty(A^^„, in the sum- 
mation. Next, we derive this bound. 

Because sinc(a;) decreases as x increases for |x| < | 

iy(A(^ > y^r^sinc (^^) . (G.15) 

Applying this inequality, 

\ 2 

A//2-1 

J2 H^(A«„,;9<'0cos(x^') 

^n=-AI/2 

, . / M/2-1 



> 1 - 

M 



J2 sine cos (nU~a^^^ - a(")nT(l + ^)) (G.16) 

. n=-M/2 / 

As we did in Appendix |El we approximate sinc(a;) ~ cos(i?x), so that we can write 
( ICTTHl) as 

\ 2 

M/2-1 
^n=-M/2 

A//2-1 

/ i \ / 

> 



( 1 - M 5Z (^) (^/^(^'' - + T=i(o )) (G-17) 

^ ^ \n=-M/2 / 

^ ^ \ n=-M/2 

X COS (7r/,(a<" - a«)nT(l + j^)) j (G.18) 
fl- — j E cos (^C^n + cos (^^("72 H , (G.19) 



90 



where 







-5a« 
1 - a(0 








^(0 




7r/,(a'''-a«)ni 



(G.20) 
(G.21) 

^ ^ (G.22) 



1 - aW- 

Equation fIG.lQp can be recognized as a Riemann-sum approximation of the following 
continuous integral 

1 \ / /"^^/^ 



1 



A//2 



2 

sm ' 



sin (d''') ) cos (^t'^t) rft) (G.23) 



Equation (1G.23P follow from the trig identity cos(x + ?/) = cos(x) cos(?/) — sin(x) sin(t/) 
and (]E.23P follows from the fact that /*^^/2 ^^'^ (^^^'''1^ j cos ^-E'^'tj c/t = 0. We can now 



use the trig identity cos(x) cos(t/) = | cos(x + y) + \ cos(x — y) to rewrite (1G.24P as 
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follows 



1 - 



M 
1 



M/2 



-M/2 



COS C<-'H cos cos E^'H dt 



1 - 

1 
1 



1 

M 
1 

M 
1 



COs2(D«) 
COs2(D(") 



M/2 



M/2 



\ COS ( (E« + &'^)t] + i COS ( (E« - C^OO 



(B(0+C'(0)fj sin (^(£W-C'W)f j 



1 \ M2 



^ .^2sinc(^(E« + C(")f j+f sinc((^«-C'«)f))' 



M 



cos^pc'^ 



X ( sine ( (E« + C"")^) + sine {{E''' - C'^'O^) 



(G.25) 



Finally, by combining flG.7p , fIG.lQP , and flG.251) and plugging back in defined expres- 
sions for C<-'\D'^'\ and (equations (10201) . (^J^, and (10221) ). we find 



(0^ 



M/2-1 
n=~M/2 

> (i-mJ — cos2(f(7'"-7'")) 
X lysine [jyafM'^ - a^'^)T (^1 + 

+ sincf (7r/,(a(')-a''))T( 1 + 



1 - a(" 
1 



1 - a('' 



1 - a(" y 27ry 



1 - a(') / 27ry 



(G.26) 
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Since cos(a;) and sinc(a;) are symmetric around x — and we have assumed 

1^(0 _ q(0| < ^ a^nj 1^(0 _ < (.^n write 

M/2-1 ^ 
n=-M/2 
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APPENDIX H 
Upper bound on Cj 



In this appendix, we derive an expression u;''' (defined in ( 15.22^ ) in terms of M . 
To do tliis we first note that, using (ID.llI) . we can upperbound a;*'' as follows 



A//2-1 Af/2-1 

,(0 A 

M2 

r=-M/2 e^r=-M/2 



^ A//2-1 Af/2-1 ^2 
r=-M/2 e^r=-A//2 ^ ^^e'*^ P 



Using our assumption a*'^ < ^ and a*'' < we can approximate -st^ 



1 + a'*' for large M. Using this approximation, we can express A^'^^^ (defined in ( ]5.17p ) 
as follows 

A« = rT-eT-a«eT-T5«, (H.3) 
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where 5''' = ^{Y'^ ~7'''} ^ j- Using the change of variables m = e — r and combining 
fHOj) and fllLSl) . it follows that 

M/2-1 A//2-1 

c:;« < ^ y y ^ (H.4) 

M-1 M/2-1 ^ 



= ^y y 

-1 M/2+m 

Af-1 A//2-1 

< — ^E E — (H.6) 



Notice that, since |e| < ^, 5''^ < and a*'' < then 



In this case 



_a«e + 5«| < |_a'''e| + |5(')| = |a«| |e| + |5«| (H.7) 



M-1 

cD''^ < ^Et (H.9) 

~ M7r2/?(0 (m - 0.5)2 

/ M-2 \ 

(H.n) 

V m=l / 

< ?<l±lB (H.12) 



M7r2/5 



using the fact that — - < —. Finally, since 



m=l 

1 1 



(H.13) 

^ l + a<'> (H.14) 
< 1 + ^ = ?^l±i, (H.15) 
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we have 



.„, ^ 2(4H-.76) 2M+1 

M7r2 2M ^ ^ 



Now if we say that M > 10, we know that < §, so that 

2(4 + 7rV6) 21 1.2011 



^(0 < ^ ' / ^ _ < _ . (H.17) 
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